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The ratio of the il^' over the J/V^ production cross section in the dielectron channel has 
been measured in ^/s — 200 GeV p -\- p colhsions with the PHENIX detector at RHIC. The 
analysis is based on fitting of the dielectron invariant mass spectra in the area around the 
J/'0 and ^1^' signals in order to extract a ratio ^|:' over J/'0 of 0.019it0.005(stat)±0.002(sys) 
and a fractional feed-down contribution to J /^|: from ij;' of 8.6 ± 2.5%. 



I. INTRODUCTION 

J /il) production in hadronic collisions can occur in part through the production of higher excited 
resonances, ip' ^ Xc which decay into the ground state. For J/^l) it is known experimentally that 
about 40% of hadroproduction rate comes from feed-down of higher mass resonances [1], [2], [3|. 
Measurement of excited charmonium states is of key importance to account for the feed-down 
contributions. Thus, this is an aditional test for different production mechanisms of charmonium. 



II. MEASUREMENT METHOD 

The number of events in the J/ip peak, Nj/^^ is given by: 

Nj/^ = ct{J/^) • B{J/^ ^ /+r ) • £ • e, (1) 

i.e. the product of the cross section {a{J/ilj))^ the branching ratio into dilepton pairs (S), the 
integrated luminosity C and the total reconstruction and trigger efficiencies, as well as acceptance 
(e). The luminosity is identical for all charmonium states, and therefore cancels in the ratios. The 
ratio of the ?/^' and J/?/^ cross sections in the e+e~ channel, R^'{e) is equal to: 

^^'(ej — (2j 

^ • o- Mj/^ e' 

where e denotes the leptonic decay channel, is the J/il){il)') production cross section and 

B{B') is the branching ratio for the e+e~ decay into the J/ipi^l)') meson. 

The ratio jfjj-^ or 'raw -jp^ ratio' is defined from the fit of the J/ip and ip' signals and must be 

corrected by the efficiency ratio, -p-, where e is defined above and e' is the corresponding efficiency 
for detection of mesons. The efficiencies are evaluated by simulation as described inllVl 



III. EXPERIMENTAL SETUP 

Data collected during the 2006 RHIC run with the PHENIX central arms spectrometers which 
cover 1 7^ I < 0.35 in pseudorapidity and 2 X7r/2 in azimuth were used for this letter. Electrons and 
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positrons are reconstructed in the central arms [4] using Drift Chambers (DC) and Pad Chambers 
(PC), located outside an axial magnetic field, and are identified by hits in the Ring Imaging 
Cherenkov detector (RICH) and by matching the momentum with the energy measured in an 
Electromagnetic Calorimeter (EmCal). 

This analysis is based on collision events triggered by minimum activity in the Beam Beam 
Counters at 3.0 < |7^| < 3.9 (minimum bias trigger), and associated hits between the RICH and a 
400 MeV cluster in the EmCal. The data sample cosisted of 88 billion minimum bias events. 



IV. ANALYSIS AND DISCUSSION 



Two sets of electron identification cuts were used to count J/il) and ?/^^ A first set, named 
'tight' with a match between energy and momentum satisfying {E/p— 1) > -3 standard deviations 
(cr), interaction- vertex cut of ±30 cm. The second set, named 'loose' with energy- momentum 
requirement (E/p-1) > -4cr, interaction- vertex cut of ±30 cm and a minimum number of DC hits 
used to reconstruct each track of the electron-positron pair. 

The number of J / candidates is obtained by counting the unlike sign pairs in the mass window 
(2.88-3.32) GeV/c^ and the i/j' candidates in the mass window (3.48-3.92) GeV/c^ The combina- 
torial background subtraction is performed as the sum of the like sign pairs. The subtracted plot 
contains charmonium resonances and physical background made of correlated (e+e~) pairs from 
Drell Yan and leptonic decays of charm and bottom. From PYTHIA the (e+e~) invariant mass 
distribution was obtained for all sources of physical background under the J/i/j and ip^ as shown 
in Figure 1. The resulting sum of all the three components was fit with a power law. 

To count the number of J/ips and t/j's different fit procedures were applied for each one of the 
electron identification cuts. Figure 2 shows a double gaussian shape for the J/ip peak, a single 
gaussian for i/j' . The background was fit with a power law, the J/i/j-i/j^ mass difference was fixed with 
the Particle Data Group value (0.589 GeV/c^) and both J/tJj and tp' widths varied independently. 
For this fit, loose electron identification cuts set was applied. 
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FIG. 1: e^e invariant mass distribution for all 
sources of physical background under J/ip and ip' . 
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FIG. 2: Invariant mass distribution e e~ in the 
region of the J/ and ip' peaks for triggered dielec- 
tron events. 



Acceptance and reconstruction efficiency were calculated by generating J/i/j and i/j^ with 
PYTHIA event generator and checking the response of PHENIX detector with a GEANT based 
Monte Carlo. In order to generate J/V^, all J/tJj decays were turned off except the dielectron one. 
To account for direct and indirect production (through x and i/j' decays) the subprocesses reported 
on the left of Table I were used. In order to generate i/j' particle id= 100443 was requested, instead 
of 443 and all ip' decay modes were turned off except the dielectron one plus the following on the 
right of Table I. 

The average transverse momentum is about 1.7 GeV for J/ip and 1.9 GeV for tp' as shown in 
Figures 3 and 4. The electron average transverse momentum is about 1.5 GeV for electrons coming 
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TABLE I: J/^ and ^|:' subprocesses in PYTHIA event generator. 




FIG. 3: Transverse momentum distribution for all FIG. 4: Transverse momentum distribution for all 
J/ijj with rapidity —0.5 < y < 0.5 generated with with rapidity —0.5 < y < 0.5 generated with 

PYTHIA event generator. PYTHIA event generator. 



from J/tp decays and 1.8 GeV for electrons coming 

J/i/j and i/j^ with rapidity —0.5 < y < 0.5 were selected and the two electrons from their 
decay were put through the PHENIX simulation chain (PISA and PISAtoDST). After that, the 
reconstructed dielectron invariant mass was determined and the electron identification cuts used 
in the data analysis were applied to evaluate the acceptance and reconstruction efficiency. The 
acceptance x reconstruction efficiency for J ftp and tp^ are shown in Figures 5 and 6 respectively, 
both for loose electron identification cuts. The integrated acceptance x reconstruction efficiency 
is 2.48 for ip' and 2.41 for J/ip. 

Integrating over the momentum distribution predicted by PYTHIA we find that the ratio of the 
ip^ acceptance x reconstruction efficiency to the J/ip is e'/e = 1.03 ± 0.01. 



V. SUMMARY 
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FIG. 5: The acceptance x reconstruction efficiency 
for J/ip for loose electron id cuts. 



FIG. 6: The acceptance x reconstruction efficiency 
for ip^ for loose electron id cuts. 



The ijj' over J / ip ratio is quoted as obtained 
with loose electron identification cuts and the 
fit double gaussian (for J/V^), single gaussian 
(for ip^) and power law (for the physical back- 
ground) -e'/e (acceptance x reconstruction ef- 
ficiency): 0.019 ± 0.005. The systematic error 
derived by taking the maximum- minimum in- 
terval for the ensemble of the various fit results 
and divided by a/T2 is ±0.002. Figure 9 shows 
the ratio results of the experiments listed on 
Table II, and also the result quoted in this let- 
ter. Based on the results here presented, the 
fractional feed-down contribution to J/ip from 
ip' is 8.6 ± 2.5% in good agreement with the 
calculation provided by [s], which is 8 ± 2%. 
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FIG. 7: ip' over J/ip ratio as a function of the en- 
ergy in the center of mass frame measured in various 
experiments. 



TABLE II: over J/ip ratio for different energy regimes and collision species. 
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